Bragg resonances for tunneling between edges of a 2D Quantum 

Hall system 

M. B. Hastings and L. S. Levitov 
Physics Department, Massachusetts Institute of Technology 

Abstract 

A theory is presented for tunneling between compressible regions on the sides 
of a narrow incompressible Quantum Hall strip. Assuming that electron inter- 
actions lead to formation of a Wigner crystal on the edges of the compressible 
regions, we consider the situation when the non-conservation of electron mo- 
mentum required for transport is provided, in the absence of disorder, by 
umklapp scattering on the crystal. The momentum given to the crystal is 
quantized due to the Bragg condition, which leads to resonances in tunneling 
conductivity as a function of the incompressible strip width, similar to those 
reported recently by N. Zhitenev, M. Brodsky, R. Ashoori, and M. Melloch. 
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The structure of the edge of a Quantum Hall system is central for understanding the 
QH physics [|T|. Theories of the QH edge predict that the edge states form an interesting 
strongly interacting one-dimensional system. Different aspects of the edge dynamics are 
described by the chiral Luttinger liquid theories , and by theories taking into account the 
long-range character of electron interaction exchange [Q and correlation effects [||. One 
expects that the QH edge has a rich structure. To probe it, one can study electric transport 
across the edge, e.g., by tunneling into the edge ||, or by measuring an ac response. The 
latter was realized in a recent ac capacitance measurement 0], in which a capacitance probe 
attached above the 2D gas was used to accumulate charge under the probe and to form a 
sharp gradient of electron density near the probe edge. When the densities in the 2D gas 
and under the probe are greater and less than an integer, correspondingly, an incompressible 
strip is formed near the probe edge. Since electric field near the edge is strong, it was possible 
to create a very narrow incompressible strip of width estimated as 8-10 magnetic lengths. 

Measuring the frequency dependent ac capacitance gives tunneling conductivity across 
the strip. As a function of the potential V g on the probe, the conductivity displays res- 
onances, uniformly spaced in V g 0. The resonances were studied at different magnetic 
fields and temperatures. The resonances shift as the field changes, however, the shift is 
much weaker than corresponding change of filling fraction, i.e., the resonances are more 
sensitive to V g than to the field. It was found that the resonances broaden with increasing 
temperature, and that the peak values of conductivity increase roughly as T a , with a > 1. 

It was demonstrated that the peaks are not caused by resonance tunneling through 
impurities: the number of peaks is about the same for the probes of different sizes, and the 
conductivity at the peaks scales as the probe size. The high degree of the peaks ordering, 
as well as the temperature dependence, inconsistent with resonance tunneling, suggests 
looking for a different mechanism. Since the strip is narrow, it may be possible to tunnel 
across it without scattering off of impurities. Momentum conservation then requires that 
the tunneling electron give up momentum to other electrons. If the other electrons (or holes 
in the region of less than integer filling factor) form a Wigner crystal [^,10[| , momentum can 



only be given up at multiples of the inverse lattice spacing. This effect would lead to Bragg 
resonances in tunneling conductivity when these multiples coincide with the momentum 
required to tunnel across the strip. 

In all theories of the QH edge predicting structure of a few magnetic lengths wide |3] [| . 
there is room for a ID Wigner crystal. For example, the theory of edge reconstruction || 
leads to an edge composed of Luttinger liquids, one chiral, and several non-chiral, each of 
the latter representing nothing but a ID crystal. One expects no long-range order in a ID 
system, so it is equally correct to view a ID crystal as a correlated fluid. The correlation 
length, however, varying as the inverse of temperature, becomes large at low temperatures. 
This leads to finite widths of the peaks in the structure factor, scaling with temperature. 

We will discuss a Wigner crystal with varying density, formed near the incompressible 
edge, and study tunneling into such a crystal. For the edge structure we use the long-range 
screening picture ||, however, since only one row of the crystal will be important, our 
predictions will hold for any reconstructed edge. We consider Bragg resonances in tunneling 
conductivity, and calculate the dynamical Debye- Waller factor due to the crystal zero-point 
and thermal motion, giving the temperature dependence of the peaks. The peak intensities, 
thermal broadening, and shifts caused by varying temperature, are found to be consistent 
with the experimental observations Jt|. 

Wigner crystal at the edge: In magnetic field, electrons form a crystal at sufficiently low 
densities and temperatures, where potential energy overcomes kinetic energy ||11|| . The 
critical density v was estimated to be slightly below 1/5 [0. At higher densities a QH fluid 
will form. 

On the sides of an incompressible strip, formed by one or several fully filled Landau levels 
with integer density, the electron density is non-integer and varies smoothly on a scale, set by 
long-range electric forces, of many magnetic lengths [[J . Thus the density of excess electrons 
and holes right near the strip must be quite small. Within the two small density regions 
there are several crystalline rows of electrons and holes, respectively (see Fig. |T].b). The 
number of the rows is determined by the external electric field E. In the continuum limit 
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|3j, the electron density measured from the integer within the strip goes to zero at some 
point with a square-root singularity. The discrete crystal will have its front row located near 
that point. The rows separation can be determined by balancing forces, eE from confining 
potential against the force e 2 /er 2 from neighboring particles, which, up to a factor ~ 1, gives 



the spacing in the front row: a = ye 2 /eE. In the experiment |7[], the e- and h-crystals both 
have small density at the strip edge, and increase away from it. We estimate that a is of 
the order of 2-3 magnetic lengths for the front row of each crystal. In the continuum limit 
||, the particle spacing goes as the —\ power of the distance from the edge. Similarly, in 
numerical simulations |§ , the particle spacing varies very slowly with the distance from the 
edge. However, the estimated density of the front row is close to the Wigner crystal critical 
density |12| , and so we expect that beyond a few crystal rows a QH fluid is formed (Fig. [j].b 



shows e- and h-rows obtained by a simulation similar to 0). 

At constant density, a Wigner crystal forms a perfect triangular lattice. In the presence 
of the field forming the strip, the density of the crystal will increase along the field, however, 
because of the long-range interaction, the crystal forms well-ordered rows, even at the edge 
||. The largest simulations of a crystal in an external field that we know of were done 
by Pieranski using 1 /r 3 -interacting dipoles [TB| (see Fig. |TJ.a for a picture of the resulting 
crystal). These experiments also revealed the remarkable idea of the conformal crystal ||13|| . 
Since the electron interaction 1/r is longer-range than 1/r 3 , the front rows in the 1/r crystal 
are ordered even better than in the 1/r 3 crystal (see Fig. 2 in ||]). 

We expect that both electrons and holes will be surrounded by spin textures, topological 



excitations of the local spin ||14|| . We will ignore any effect this has on the crystal, as well as 
any question of how the spin is transferred across the strip. It may be that spin is transferred 
independently of charge, by spin waves. 

Bragg reflection Since the tunneling amplitude falls exponentially with distance, the main 
contribution to the tunneling rate arises from transitions between the front e- and h-rows 
(Fig. [TJ.b). Therefore, we will focus attention on only one row on each side of the strip. We 
consider a strip of width w going in the x— direction, with e- and h-crystal row spacing a. 



In Landau gauge, electron states can be written as 

iP{x,y) = e i v*f n (y+pcK/eB) (1) 

where /„ is the harmonic oscillator n— th wave function, p labels the x-momentum. 

The tunneling amplitude is given by the overlap of states (]!]) on opposite sides of the 
strip. In the absence of momentum non-conserving scattering, the amplitude will vanish 
because shifting y by w across the strip is equivalent to a change of p, but the states with 
different p's are orthogonal. In the presence of a crystal, electron momentum may change by 
a multiple of the Bragg vector, Sp = 27m/a. (For simplicity, let us assume that the crystal 
period is the same on both sides of the strip.) On the other hand, since according to ([!]) the 
change in momentum is coupled to a y— shift, the overlap of the states will be maximal at 
5pch/eB = w. The two conditions put together give that tunneling must be maximal at 

™n=^n, (2) 

where n corresponds to the Bragg resonance order, and $o = hc/e is the flux quantum. 

From (fj), the peak positions are dependent upon the strip width w, and thus are very 
sensitive to the gate voltage V g . Increasing V g leads to squeezing of the strip: roughly, w 
is inversely proportional to V g . For this reason, we expect that the dependence on V g of 
the B— field needed to observe a given resonance is stronger than the dependence on V g of 
the B— field needed to maintain a given filling fraction under the gate. This is seen in the 
experiment. 

To estimate the order of a resonance one takes the product wa/l H , where Ih is magnetic 
length. In the experiment, Ih — 10nm, w is about 5 — 8 Ih, and a is estimated above 
to be 2 — 3 Ih- This gives the order of resonance between 10 and 25, consistent with the 
experiment. 

Going from a 2D crystal to single row: To estimate the response of the crystal to injecting 
an electron, we assume that all relaxation occurs within one row, and treat it as a one- 
dimensional system. The other electrons will be taken into account only to the extent that 
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they produce a confining potential V(y) for that row. This approximation underestimates 
the crystal stiffness, enhanced by the presence of other rows, and thus overestimates the 
peak widths. However, we believe the discrepancy is not large, because for a crystal with 
long-range interaction the shear modulus is much smaller than the compressibility, and also 
the particles in the front row, being more widely spaced, are more loose than in other rows. 
Thus, charge spreads most easily within the first row. 
With this, one writes the Hamiltonian for one row as 



where r n = x n i + y n j are electron positions, and the interaction U(r) = e 2 /e\r\ — e 2 /e(r 2 + 



In a magnetic field, kinetic energy is very small compared to potential energy, and so 
the fluctuations about mean positions in the crystal will be small. We assume that different 
Landau levels are not mixed by the interaction U, and then treat the displacements x n 
and y n as conjugate variables with the commutation relation [x n ,y n ] = i$ /B. Expanding 
the Hamiltonian, and keeping only quadratic terms, one obtains a harmonic chain problem. 
Diagonalizing the Hamiltonian, one gets elementary excitations. Their spectrum at small 



qa <C 1 is acoustic: uj(q) = c\q\, where c = ayre/m*. (Here k = J2n n2 U" =a ni m * = 

m(l + muJc/V" =0 ).) At small q, the displacements are along the row. 

Evaluating the tunneling rate: The tunneling is described by transferring an electron across 
the strip, from the e-row to the h-row. Injecting an electron into the row is accompanied 
by motion of other electrons that give space to the new electron. This makes the tunneling 
a collective process, and requires considering relaxation in the row. Since at high magnetic 
field the effective mass of electrons in the row is large, m* ^> m, the relaxation will be slow. 
This simplifies analytic treatment, making it possible to consider only the effect of acoustic 
modes, qa <C 1. In order to evaluate the tunneling rate we assume that the transfer of an 
electron across the strip is a much faster process than charge relaxation within the row, 
and that the characteristic time scales are much larger than a/c. Such an assumption is 




rf 2 ) 1 / 2 accounts for the screening by the gate. 
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self-consistent at low T < hu c , which is within experimental range (see similar discussion 
in 

We evaluate the tunneling rate semiclassically, by taking the saddle-point of the action 
for the dynamics in imaginary time corresponding to injection and removal of an electron 
in the rows. The path in imaginary time must be closed ("a bounce"), evolving the system 
to the initial state ||15|| . The tunneling rate is thus given by the sum over all injection and 
removal events, (ri,ti), (r2,t 2 ), with appropriate phase factors: 

G(w) = W V[u^\ u W] e ^(n-r 2 ) e -s[f-s^ ^ (3) 

1,2 •* 

where 1 and 2 label different positions in the chain, w = 2^Bw/^q, the labels (e) and 
(h) mark the electron and hole rows, and ■u (e,/l ) is the displacement field in the rows. It is 
convenient to rewrite the tunneling conductivity G(w) as a convolution, 

G(w) = J F {e \w - s)F {h \s)ds , (4) 

where F^ e ' and F^ are the structure factors of the rows: 

F(w) = E / ^Me^ (ri - r2) e" 512 . (5) 

1,2 

The usefulness of using F(w) is that the e- and h-rows may have different densities or elastic 
moduli, which will make spacings and widths of peaks of the conductivity G(w) irregular, 
while the structure factor of each row still will be a simple function. 

Let us calculate the structure factor of one row. Because of what has been said above we 
consider only the long-wavelength displacements u(x,t). Choosing the units so that a = 1 
and c = 1, the action for one row is 

S = J J (^u 2 + ^u 2 x + X x u - Xu^j dx dt , (6) 

where p = m*/a, and the Lagrange multiplier X(x,t) is a multivalued function changing by 
±2iri going around the points (x%,ti), (x2,t 2 ), respectively. The term with A is introduced 
to describe injecting an electron at (xi,ti), and removing it at (^2,^2)- Qualitatively, the 
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particles are moving as shown in Fig. |T].c. The equation of motion then is V 2 w = 0, with 
boundary conditions: 

a) The displacement u(x,t) changes by ±1 going around the injection and removal points; 

b) As a function of time, u(x, t) is periodic, with period (3 = 1/T. 
The solution is given by 

T 1 i sinh7rT(z - z x ) 

u (z) = Im — In , {- , 7 

z7T smh Til (z — z%) 



where z = x + it. Evaluating the action S[uo{x, t)} one gets ^ In | sinhf — ^ 
is equal to 



which 



5 = -ln-(smh _ + sm — ) , (8) 

where X\2 = X\ — X2, t\i =t\ — ti. 

To average the phase factor in we write ri^ in terms of the displacement u(xi,2) from 
the mean position explicitly showing the particle number in the row: 

F(w) = V] ^_-jn^ e iii(n+«(xi,ti)-u(a:2,t2))^ e -S , o ^ (g) 
H2=n 

where the sign is due to exchange of n + 1 fermions, and the Debye- Waller factor (...) 
describes dephasing due to the zero-point and thermal fluctuations. The gaussian average 
in (H) is performed in a standard fashion, 

(...) = e iu>n e -\w 2 {(u(x 1 -T/2)-u{x 2 ,T/2)) 2 ) (10) 

where the second exponent equals —(wh/pac) 2 So. Taking the extremum of (|J) in the time 
separation t 12 gives optimal t i2 = (3/2. Finally, restoring c and a, the resulting structure 
factor is 

where T = 2naT/hc, a(w) = pa 2 c/2nh + w 2 h/2irpc. The structure factor (pT|) has peaks at 
w m = (2tt /a)m, in agreement with (Q). The width of the m— th peak scaled in a -1 is given 
by a 1 / 2 (w m ) T, which equals 



[7 + m 2 /^} 1/2 f , 7 = pa 2 c/27ch . (12) 

The peak width increases with m and is proportional to temperature. 

To compare with the experiment [J], we plot the ac capacitance 
C(w) = (7/(1 + {Cuo/G(w)) 2 ) with G given by (g), (0)' and ^ bein S a parameter (Fig. f2|). 
The peaks appear at the transition between plateaus, and have temperature dependence 
consistent with fTj. 

Qualitative discussion and conclusion The distinction between the Bragg resonances and 
other explainations of the observed peaks can be made by studying thermal shifts of the 
peaks. In the Bragg condition @ both the strip width w and the crystal period a depend on 
T, which results in the shifts of the same sign for all the peaks, and of magnitude proportional 
to the order of the peak. The sign of the shifts is determined by the competition of two 
factors. The thermal expansion of the crystal will lead to a increasing, and w decreasing 
with temperature. According to (0), the first effect is more important, which leads to a shift 
consistent with experiment. 

The effect of disorder is difficult to estimate, but in the setup [0] the sensitivity to disorder 
should be lower than in other existing devices: there is a stronger electric field, and thus a 
more narrow strip; both sides of the strip are free, and so the disorder potential will mainly 
result in the strip bending with lesser effect on the width. 

We may conclude by listing several qualitative features of the experiment correctly de- 
scribed by the above theory. The number of peaks and their positions in the magnetic 
field-gate voltage plane are a result of the Bragg condition. The change in the height of 
peaks is a result of the temperature dependence of the action (§). The shifts in peak position 
are explained by a thermal expansion of the crystal. 
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FIGURES 

FIG. 1. a) A large collection of magnetic dipoles which repelled each other was confined 



between two glass plates, and the glass was tilted; the in-page gravity component is down [13] 
(courtesy of P. Pieranski and F. Rothen). b) Incompressible strip with few electron and hole rows 
on the sides is shown. Beyond the rows are regions of QH fluid. The coordinate x is along the strip, 
and y is normal to it. c) The imaginary time dynamics picture. The world lines of the tunneling 
particle and of other particles are shown. 



FIG. 2. The ac capacitance C(w) = (7/(1 + (Cuj/G(w)) 2 ), with G given by (|) and (0), 
plotted as function of the strip width w = 2ttwB/^q at several temperatures, offset for clarity. 
Parameters used: Cuj = 0.05, = 1, a" = 1.68, 7 = 2. Non-equal e- and h-periods are chosen 
to demonstrate the effect of possible period mismatch. 
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